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Abstract Diffusive shock acceleration is the theory of particle acceleration through 
multiple shock crossings. In order for this process to proceed at a rate that can be 
reconciled with observations of high-energy electrons in the vicinity of the shock, 
and for cosmic rays protons to be accelerated to energies up to observed galactic 
values, significant magnetic field amplification is required. In this review we will 
discuss various theories on how magnetic field amplification can proceed in the 
presence of a cosmic ray population. On both short and long scales, cosmic ray 
streaming can induce instabilities that act to amplify the magnetic field. Develop- 
ments in this area that have occurred over the past decade are the main focus of 
Oh. this paper. 
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1 Introduction 

^ , The theo ry of diffusive shock acceleration (DSA) originates from the idea originally 

posed bv iFermil ( 1949(1 that cosmic rays are scattered by waves to be isotropic in 



their local frame. Scattering was proposed to happen in clouds where a turbulent 
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magnetic field would isotropise the cosmic rays. Since on average head-on colli- 
sions are more frequent, a net acceleration occurs that naturally creates a power- 
law. Diffusive shock acceleration, proposed by various authors in the late seven- 
ties (JAxford et alJll977t lKrvmskiilll977l lBe!illl978alll978ri iBlandford and Ostrikerl 



Il978f) is the discovery that this acceleration proceeds much faster in the vicinity 
of a shock. When crossing the shock, the first collision is always head on, thus 
allowing the acceleration to proceed at a significantly faster rate, making it an at- 
tractive mechanism to accelerate cosmic rays to high energies. The faster the shock 
velocity, the larger the energy gains upon transition of the shock. The theory pre- 
dicts a powerlaw with a spectral index that matches quite closely the observed 
powerlaw of cosmic rays arriving on Earth. 

The first argument that pointed to supernova remnants ( SNRs) as the prime 
sourc es of Galactic cosmic rays was based on the energy budget (jBaade and Zwickvl 
ll934T) . The energy required to replenish the cosmic rays against their losses from 
the Galaxy amounts to about 10% of that available in SNRs. However, it was not 
until the late seventies that a theory on how to transfer kinetic energy efficiently 
into the acceleration of cosmic rays had been developed. As shown in this review, 
even though the basic test particle theory is simple, the intrinsic nonlinearity of 
this process makes it a difficult one to grasp in full and more work needs to be 
done to understand the process from start to finish. 

The strong support that DSA can work as the major mechanism to accelerate 
cosmic rays came from observations of the thin X-ray rims at the blast waves of 
SNRs. The magnetic field required to confine the cosmic rays in the vicinity of the 
shock is much higher than the mean interstella r magnetic field ( Achterberg et al.l 



11994 IVink and Laming! l2003t IVolk et al.l 120051) Originally it was proposed that 
the magnetic field could be amplifie d resonantly ( Lerchefl967tlKulsrud and Pearcd 



Il969l: IWentzelll974l: ISkillindll975bl) . but it is not clear that the magnetic field will 
be able to grow beyond 5B/Bq ps 1 when the resonance condition is lost, as would 
be required in order to explain the observations. In the past decade, a number 
of theories has arisen that could potentially explain amplification of the magnetic 
field to values corresponding to those observed in supernova remnants. 

However, there is no reason to believe that only SNR shocks accelerate parti- 
cles. Shocks are abundant in the universe on all scales, and how and when they 
become efficient accelerators is an active area of current research. Locally, cos- 
mic rays are accelerated in heliospheric plasmas. On larger scales, evidence ex- 
ists for_accelerated elect rons i n the lobes of radio galaxies, as suggested early on 
by IBlandford and Reed (11974) and backed up by more detailed observations an d 
modelling fe.g. lCarilli et al.lll99lt ICroston et ai]l2009t iBlundell and FabiarJlJOllh . 



Furthermore clusters of galaxies have been obs erved to contain nonthermal radio 
emission suggesting active particle acceleration iBagchi et al.l (J2006I ) : iFerrari et al.l 
(|2008h . 

Support for DSA on much larger scales has only recently been discovered with 
detailed and s pectral obs ervations in radio wavelengths of Megaparsec scale shocks 
( van Weeren et al.ll2010l) . The spectral hardening downstream from the supposed 



shock front in these systems indicates that active acceleration proceeds at the 
shock front itself. The shocks are believed to be the result of mergers of clusters of 
galaxies. The magnetic field deduced from the radio observations indicate a field 
strength of the order of /iG, which is much higher than would normally be expected 
in the intra-cluster medium ( Briiggen et alj|201ll) . This is an additional indication 
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the the process of DSA and magnetic field amplification is intrinsically linked and 
occurs at shocks on all scales. Because of the scales of these shocks potentially 
protons can be accelerated to energies of 10 19 eV, although no direct evidence has 
been found yet. Magnetic field amplification through streaming cosmic rays has 



also been suggested as a source for primordial magnetic fields ( Miniati and Belli 
20111) . 

On the Galactic scale, in addi tion to SNRs, shocks around superbubbles have 
been discussed as accelera t ors (e.g. Bvkov and Fleishman 1992: Bvko v and Topt ygin 



2001 



; IParizot et all 120041: iBinns et al.ll2008t iButtl l2009t iFerrand and Marcowithl 
20101) . Superbubbles are formed as a result of cumulative outflows from an assem- 
bly of massive stars, possibly enforced by the supernova explosions themselves. 
The outer region thus formed is a shock of large proportions that could poten- 
tially be a cosmic ray accelerator, or various multiple shocks can consecutively act 
to accelerate particles, which may modify the resulting spectral index. 

Although in reality many deviations arise due to e.g. nonlinear modification of 
the shock structure, magnetic field obliquity, geometric effects, time-dependence, 
and magnetic field amplification, the basic theory still holds. However, in order 
to be able to compare the theory with observations, all of the complicating fac- 
tors need to be taken into account. In this review, we focus on the advances in 
theory that have been made over the past decade. More specifically, we will dis- 
cuss how the amplified magnetic field required for efficient acceleration is a direct 
result of DSA and the presence of co smic rays. Details of the other processes 



can be found in earlier review articles (|Drurvlll983l : iBlandford and Eichlerlll987l : 



iMalkov and O'C DrurvlboOll [ffillas 2005) 



In Section [2] we will briefly review the original theory of diffusive shock accel- 
eration and demonstrate how the process naturally results in a power law cosmic 
ray spectrum. In Sections [3H5] we will discuss a number of theories that couple dif- 
fusive shock acceleration to magnetic field amplification. We will mainly focus on 
the most recent theories. In Section [6] we will discuss possible deviations from the 
source spectrum as a result of shock obliquity, nonlinearity, and time-dependence. 
We refer to other chapters in this book for a treatments of DSA in relativistic 
shocks, and to the observational insights and developments. 



2 Diffusive shock acceleration 

The powerlaw in energy that results from diffusive shock acceleration can be under- 
stood in different ways, highlighted by different authors in its period of discovery. 
A crucial point is that the cosmic rays isotropise on either side of the shock due 
to small-angle scattering off magnetic field fluctuations. The faster the isotropisa- 
tion, the faster the particle can recross the shock. Every time the shock is crossed, 
a net energy gain is received by the particle crossing the shock. Although the 
acceleration efficiency depends on the effective scattering efficiency, the resulting 
spectrum is independent of the diffusion coefficient. In many cases the most effi- 
cient scattering rate of Bohm diffusion, where the mean free path is of the order 
of the gyroradius, is used in order to generate cosmic rays with the high energies 
that are observed. 

An intuitive way of approaching the problem is by evaluating the number of 
particles that is located at the shock versus the number of particles that escape 



K. M. Schure et al. 



downstream. Only particles that do not escape qu alify for the next round of acceler- 
ation. This is the approach originally described by Bell ( 1978a|) and described as the 
microscopic approach in the review by Drur vl (|l983|). T h e macroscopic approach, 
originating from iKrvmskiil (|l977l) : lAxford et all (|l977t ); iBlandford and Ostrikerl 
(jl978| l. derives the acceleration and resulting powerlaw from the distribution func- 
tion, requiring continuity at the shock. Below we will briefly summarise both meth- 
ods and we refer to the original papers or iDrurvl (1983) for a more extensive 
treatment. 



The distribution of relativistic particles can be described by the Vlasov-Fokker- 
Planck equation: 

df , / \df ldu df d ( . ,df\ n „. 

The distribution can be separated into an isotropic part, and anisotropic parts 
to arbitrarily high order: 

/ = /o + ^ + ^ + ... (2) 

v v 

Mostly, the diffusion approximation is used, in which the first order anisotropy 
(/i = fiv.i/v) is used and eliminated, to arrive at a distribution that depends on 
the isotropic cosmic ray density (/o) alone. 



^-+u^- + -^- - -— p^- = (3) 

dt dz 3 dz 3 dz dp 
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-"/i (4) 



dfo . df d ( c 2 df \ ldu df 



dz 

(5) 
dt " dz dz \3v dz ) 3 dz 1 " !> 

where we used k = c /(3u). 

From Eq. \5\ far upstream the steady state solution (dt = and du/dz = 0) 
implies that /o should have the form 

/ = /o(p)e~(^) z (6) 

to have a bound solution, where we use that in the shock frame u = — u\. Or 
alternatively, 

§/i = ui/o. (7) 

Downstream the steady state solution gives 

/ = MP)- (8) 

At the shock the solutions have to connect, giving the boundary condition: 

c . . Ui dfo u 2 df . . 

3 /l + T p ^ = T p ^- (9) 
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Using the boundary condition for the far upstream we can replace c/i/3 with mi/o, 
giving: 

, . /mi u 2 \ dfo , , 

which results in requiring that the cosmic ray density follows a powerlaw distri- 
bution: 

/o«P~ 9 (11) 

with 

« = 3tu/(tu-ua)=3r/(r-l), (12) 

where r = U1/U2 represents the compression ratio at the shock. This powerlaw is 
valid in the test particle approach, for a planar shock, where the magnetic field 
i s parallel to the shock normal. More details can be found in the original papers 
(|Krvmskiilll977l lAxford et ai]|l977t iBlandford and Ostrikerlll978l). 



An alternative approach was used by Bell ( 1978a) to derive the powerlaw distri- 
bution of shock accelerated particles. It is based on the microscopic physics and is 
helpful to get insight in how the powerlaw may change depending on escape prob- 
ability and probability of crossing the shock, which will be useful in understanding 
the physics of the later sections. 

The flux of particles downstream is just the number of particles n that are 
advected with the downstream flow velocity u 2 : nu 2 . The number of particles 
crossing the shock front per unit time from upstream to downstream in case of an 
isotropic distribution is half the number of particles moving towards the shock, 
and their average velocity over angle is again half of the shock velocity, giving 
for the flux nc/4. The fraction of particles not returning to the shock is therefore 
nu2/(nc/4) = 4u 2 /c. The probability of recrossing the shock can be high: P re t = 
1 - Au 2 /c. 

The energy gain of a particle crossing the shock from upstream to downstream 
can be calculated by transforming the momentum of the particle in the upstream 
to the downstream frame: p = p{\ + (ui — u 2 ) cos 9/c) such that the average change 
in momentum is 1p(u\ — u 2 )/(3c). The energy gain from downstream to upstream 
is exactly the same, as u\ and u 2 are interchanged and the angle of integration 
runs to the opposite side, yielding an extra — 1. Thus the gain of momentum 
after a complete cycle is Ap = 4p(«i — u 2 )/(3c). After k cycles, the number of 
particles has decreased as n = no(l -~4u 2 /c) and the momentum has increased as 
V = Po(l + 4(ui — u 2 )/(3c)) . The number of particles as a function of momentum 
can be found to be 

ln(n/no) _ k\n(l — 4u 2 /c) —4u 2 /c _ 3 



ln(p/p ) feln(l + 4(Mi-w 2 )/(3c) 4(wi - u 2 )/(3c) r-V 

where r again is the compression ratio r = u\/u 2 , such that: 

-3/(r-l) 



— = I JL 

no \po 



(13) 



(14) 



and the differential energy spectrum is: 

ndpo:p~ (r+2)/( - r - 1) dp. (15) 
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In terms of the distribution function we arrive at the same answer as from the 
macroscopic approach, since: 



n 



-3r/(r-l) 



A = 5? KP ' (16) 

For a more detailed treatme nt we refer to the original papers and earlier reviews 
(JBeHlll978alll978bl lDrurvlll983|) . 



3 Magnetic field amplification: resonance regime 

In the theory of diffusive shock acceleration cosmic rays are accelerated by crossing 
and recrossing a shock, as shown in the previous section. On each cycle of crossing 
recrossing between upstream and downstream the cosmic ray energy increases by 
a small fraction ~ u s /c where u s is the shock velocity. For acceleration to PeV 
energies a cosmic ray (CR) has to cross the shock ~ 10c/u s times. A shock propa- 
gating into a purely uniform magnetic field cannot accelerate CR to PeV because 
charged particles pass easily through the shock and escape upstream or down- 
stream making only one pass through the shock and gain little energy. Arguably, 
it was the realisation that charged particles are not free to escape the shock envi- 
ronment that provoked the development of the theory of shock acceleration in the 
late 1970's. 

Fast and efficient CR acceleration by the Fermi mechanism requires that parti- 
cles are multiply scattered by magnetic fluctuations in the acceleration source (e.g. 
shock) . Magnetic field amplification due to the resonant cosmic-ray streaming in- 
stability was studied in the context of galactic cosmic-ray origin and propagation 
since the 1960s (see e.g. iKulsrud and Cesarskvlll97ll : IWen tzel 1974; Achterberg 
Il98ll iBerezinskii et alJll990HZweibe]|l2003(l . It was proposed bv lBelll (|l978al) as a 



source of magnetic turbulence in the test particle DSA scenario. 

CR streaming along magnetic field lines excite unstable gro wth of Alfven waves 
with wavelengths comparable with the CR La rmor radius ( Kulsrud and Pearcd 



119691: IWentzellfl973 ISkillindll97"5alll975blll975ch . The Alfven waves consist of cir 
cularly polarised distortions to the magnetic field lines. CR gyrating along the 
field lines in spatial resonance with the fluctuations are strongly scattered and 
consequently execute a random walk along a field line. The appropriate model for 
CR transport in the shock environment is diffusion instead of free propagation. 
CR cross a shock many times with a statistical probability that naturally results 
in a E~ 2 energy spectrum for cosmic rays (|Krvmskiilll977t lAxford et al.lll977t iBelll 
Il978at iBlandford and Ostrikedll978h . as shown in Section [2] 

The theory of wave excitation and CR scattering had previously been applied 
to CR propagation through the Galaxy. Skilling set out the coupled equations for 
the wave energy density / and the CR distribution function /. When applied to a 
steady state wave and CR precursor ahead of a shock the equations take the form: 

d£_d_{±a 1 df\ dl _ v_a d(^c P 4 f/3) 

s dz~ dz\Zn I dz) s dz~B$/8n dz ' [<> 

where va is the Alfven speed, p is the CR momentum, and r g is the CR Larmor 
radius. In back-of-the-envelope terms, 47rcp 4 //3 is the CR pressure. / is the ra- 
tio of the energy density in the Alfven waves (SB 2 /Att) to the energy density of 



MFA and DSA 



the unpe rturbed f ield Bq, I = 25B 2 / Bq. Further details of the equations can be 
found in ISkilline (I1975al 1975blll975dl , and a solution of the precursor equations 



Belli ( 1978a|) . The dominant physics of the interaction between 



can be found in 

CR and the Alfven waves is on the one hand that wave growth is driven by the 
CR pressure gradient and on the other hand that the CR diffusion coefficient is 
inversely proportional to the wave energy density with mean free path A given by 
A = 4r g /3J. The equation for wave evolution can be integrated to give 

I = 2^ 4 " Cp4 / /3 (18) 

v A puj 

For a characteristic interstellar magnetic field {B ~ 3/xG) the Alfven speed is 
around va ~ 10 km s -1 , so u s /va is of the order of 10 3 for the outer shock of a 
young supernova remnant (SNR). To account for the energy density of Galactic 
CR, acceleration by SNR must be efficient and the energy transfer into CR at a 
SNR shock has to be in the range 0.1 — 0.5pu 2 ( Baade and Zwickvlll934t lLongaiii 



|2010T) . In terms of Equation [18] this gives J > 1, which would correspond to 
a perturbed field SB greatly exceeding the zeroth order field Bq- Naively this 
implies a diffusion coefficient much less than r g c and a CR scattering mean free 
path much less than the Larmor radius (A <C r g ). The linear instability depends 
upon a resonance between the CR Larmor radius and the Alfven wavelength. 
This resonance is destroyed when SB approaches Bo, in which case the instability 
cannot be resonantly driven and the instability is expected to saturate at about 
SB ~ Bq- The linear equations lose validity when I ~ 1, but Equation 1181 suggests 
that instabilities driven by CR streaming may be able to amplify magnetic field far 
beyond its initial ambient value. In the next sections we will turn to instabilities 
that do not rely on the resonance condition and may continue well beyond SB /Bq ~ 
1. 



4 Magnetic field amplification: short wavelength regime 

4.1 A non- linear estimate of the amplified magnetic field 

A more basic understanding of the opportunity for magnetic field amplification 
can be derived as follows. According to equation [T7] the rate at which the wave 
energy density I grows is at v^dPcr/dz where P cr is the CR pressure. This can 
be interpreted as a force dP C r/dz pushing against magnetic fluctuations at the 
Alfven velocity x>a- Provided the CRs are coupled to the magnetic field fluctua- 
tions and the fluctuations propagate at about the Alfven speed this equation is 
approximately valid even when the spatial resonance between the wavelength and 
the Larmor radius breaks down. 

In the linear regime it makes sense to think of the magnetic fluctuations moving 
at the Alfven velocity. In the non-linear regime, the fluctuations no longer take 
the form of linear waves, but still they can be estimated to move at a velocity 
v f ~ V A ~ (magnetic pressure/density) 1 ' 2 . Energy transfer to the fluctuations 
occurs at a rate v-fdP C r/dz, and the equation for growth of the turbulent energy 
density Uf associated with the fluctuations is 

UsdU f /dz = VfdP s /dz, (19) 
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where P s is the CR pressure at the shock. Assuming that Vf = (Uf/p) 1 ' 2 and 
Uf = B /4ir as for Alfven waves, the turbulent energy density at the shock is 



:5_~I ( Pcr 

4tt ~ 4 \pv? s 



2 



pu s (20) 



If the CR acceleration efficiency is P s ~ Q.lpu 2 this leads to an estimated 
magnetic fi eld of the order of 100 uG in youn g SNR in good agreement with ob- 



servations (|Vink and L aming 2003; Bere zhko et a l. 2003; Volk et al. 2005). For a 
fixed accelera tion efficiency P cr oc pu 2 , the amplified magnetic field is proportional 



to p 1 ' 2 u s , see lBell and Lucekl (2001) for a more detailed model. A later improved 



analysis of field amplification suggests a stronger dependence on shock velocity, 
B oc p ll2 uj 2 , as discussed in Sect. 1131 (|b1-I3|2004 laxit) . 



4.2 A non-resonant instability 

Equation 1201 sidesteps the requirement for a resonance between the Larmor ra- 
dius and the instability wavelength and subsequent theoretical developments have 
shown that the resonance is not essential for the amplification of magnetic field by 
CR streaming. The same linearisation of the Vlasov equation that gives rise to the 
res onan t insta b ility also ide ntifies a non-resonant instability as first demonstrated 
bv iBelll (|2004h . iBelll (|2004f ) simplified the problem by using the MHD equations 



to treat the thermal plasma as a magnetised fluid. A kinetic Vlasov treatment 
is retained for the CR. The CR exert a force on the MHD plasma through the 
reaction on the j cr x B force, which is equal and opposite in most of the unstable 
regime, the details of which will be discussed in Sect. 14.31 Here j cr is the CR elec- 
tric current density, and the CR trajectories are calculated in the magnetic and 
electric fields calculated by the MHD model. 

The coupled Vlasov-MHD equations give a linear dispersion relation containing 
two different instabilities for the two different circular polarisations and wavenum- 
bers k parallel to the zeroth order magnetic field. The resonant instability insta- 
bility occurs when the circular polarisation is such that CR moving in the same 
direction as j C r gyrate resonantly around the zeroth order field in the same sense 
as the helical perturbations in the magnetic field, which we will refer to as the 
left-hand polarisation. A stronger non-resonant instability is found in the opposite 
(right-hand) circular polarisation. 

The dispersion relation for the non-resonant instability is plotted in Figure [1] 
for a power-law CR distribution f(p) oc p for p± < p < P2 where f(p) is non-zero 
between momenta p\ and p2- The resonant and non-resonant instabilities have 
similar growth rates when k « r^ 1 , that is when the inverse wavenumber k 
is approximately equal to the Larmor radius of CR with the lowest momentum 
p\. The lowest momentum CR are the most important because they carry most 
of the electric current. At wavelengths less than r g there are few CR in spatial 
resonance with the helices in the magnetic field and the growth rate of the resonant 
instability decreases as k increases. In contrast, the non-resonant instability does 
not depend on resonance and its growth rate instead increases with increasing 
k. The non-resonant growth rate increases with k until tension in the magnetic 
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Fig. 1 Dispersion relation for the non-resonant instability va = 6.6 X 10 ms 1 , u a /c = 1/30. 
ui is in units of u "\ /cr g and k in un its of Tg where r g is the Larmor radius of the lowest energy 
CR. Reproduced from lBelll J2004D . 



field overpowers the driving force, that is when j cr x B = B x (k x B)c/4-7r and 
kB ~ Aixj/c. 

The essential difference between the resonant and non-resonant linear instabil- 
ities can be explained as follows. In both cases the background MHD fluid motions 
are driven by the j C r x B force. j C r and B have unperturbed zeroth order compo- 
nents, j cr o an d Bo, and first order perturbed components j cr % and Bi respectively. 
The j cr x B driving force has two first order components, j cr o x Bi and j cr i x Bo- 
The resonant instability is driven by j cr i x Bo where j cr i is the perturbed CR 
current which is especially strong when the CR trajectories respond resonantly 
to a helical magnetic field with a wavelength equal to the CR Larmor radius. On 
the other hand, the non-resonant instability is driven by the other first order force 
j cr o x Bi. In the non-resonant instability only the uniform zeroth order current j cr o 
matters and there is no requirement for a resonance with the CR Larmor radius. 

Because the non-resonant instability is driven by the zeroth order CR current, 
we can derive a growth rate while omitting the first order current from the analysis. 
The response of the CR trajectories to the perturbed fields can be ignored, and 
hence there is no need to solve the Vlasov equation for CR. It is sufficient to solve 
the MHD equations for the thermal plasma with sole addition of the j cr o x Bi 
force in the MHD momentum equation. The MHD mass conservation equation can 
also be omitted since the instability is transverse and the density is unperturbed, 
pi = 0. Similarly the thermal plasma pressure is unperturbed. The first order 
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equations are then 

^^-iBoxfVxB,) ^l = Vx(u 1 xB ) (21) 

ui can be eliminated between the equations to give 

d ~W i (Bo ■ V)2 Bl = 'Jo (B ° ' V) (JCT0 X Bl) (22) 

Harmonic solutions Bi oc exp[i(kz— cot)] with k parallel to Bo in the z direction, 
gives a dispersion relation 

2 ,22 .kBojcrO , oq n 

UI — K V A = ± (2ol 

pc 

In the absence of a CR driving current j cr o = 0, this is the dispersion relation for 
an Alfven wave. The right hand side of the dispersion relation creates an instability 
provided kBoj cr o/pc > k v A and the negative sign is chosen for ± corresponding to 
the non-resonant polarisation. If k v A > \kBoj C ro/pc\, the tension in the field lines 
overpowers the driving term and the dispersion relation is that for non-growing 
Alfven waves with a modified phase velocity. The neglect from the analysis of 
the j cr i x Bo ignores the tendency of CR to follow the field lines on wavelengths 
greater than the CR Larmor radius and therefore omits the reduced instability 
for kr g < 1 that is seen in Figure [1] For wavenumbers k in the unstable range 
rg 1 < k < 4irj cr o/Boc, equivalent to r g l < k < Boj cr o/(pv A c), the instability is 
purely growing with a growth rate 

kBojcro\ (24) 

pc ) 

and a maximum growth rate, determined by the tension in the magnetic field, 

*•— Ht)" 2 ^ (25) 

The CR current j cr o in the precursor is related to the CR pressure by 

jcrO = e(u s /c)Pcr/pi\n(p 2 /pi). 

If the CR acceleration efficiency is defined by P cr = r]pu 2 s , the number of e-foldings 
for the fastest growing mode is 

_ k yM A 
K B 21n(p 2 /Pi) 

where u> g is the CR Larmor frequency, M A is the Alfven Mach number M A = 
Anpul/B 2 , k/kb is the ratio of the CR diffusion coefficient to the Bohm coefficient, 
and T pc = (c /u s )u>J is the time taken for the shock to propagate a distance equal 
to the scaleheight of the CR precursor if Bohm diffusion applies. Since M A ~ 10 
for a young SNR propagating into an interstellar magnetic field, and overall CR 
efficiencies must be 10 — 50%, the number of instability e-folding times is much 
greater than one. The term k/kb introduces an element of self-regulation since the 
number of e-foldings is greater if the magnetic fluctuations are small and diffusion 
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is greater than Bohm. This may be important at the foot of the precursor where 
the instability has not had the opportunity to grow. If scattering is weak, CR 
escape a relatively large distance upstream, initiate instable growth far ahead of 
the shock and remedy the lack of a perturbed magnetic field able to scatter the 
CR. 



4.3 Return currents and energy transfer 

In the above analysis the force on the background thermal plasma was included 
as — j cr x B. To conserve momentum the force on the background plasma is equal 
and opposite to the force on the CR. Another way of looking at this is that the 
force on the background plasma is exerted through the return current jt that the 
background thermal plasma must carry to neutralise the current carried by the 
CR, jt ~ — jcr- The difference between jt and — j C r is given by the Maxwell equation 
when the displacement current is neglected: jt = — j cr + cV x B/4-7T in which case 
jt x B = — j cr x B — B x (V x B)c/4-7r. jt and j cr very nearly cancel out where the 
instability is strong. However, their non-cancellation at small wavelengths gives rise 
to the magnetic tension that limits the instability to the range k < Boj cr o/pv A c. 
At longer wavelengths the thermal return current follows a path very close to that 
of the CR current. 

A necessary condition for CR-driven instability is that the thermal particles 
are strongly magnetised with Larmor radii much smaller than the wavelength and 
that the CR should have a Larmor radius comparable with or larger than the 
wavelength so they are not tied to magnetic field lines. If the thermal particles 
are frozen-in to the magnetic field it is not immediately obvious how one can have 
a jt x B force with a current jt unaligned with the magnetic field. The difficulty 
is resolved through the theory of cross-field drifts. The jt x B force imparts an 
acceleration to the plasma. This acceleration can be viewed as being equivalent 
to a gravitaional force. A charged particle in a gravitational field executes a cross 
field drift, and similarly in this case a cross field drift arises through the Lorentz 
force, which produces the current density jt. 

Another conundrum is how energy is extracted from the CR current to drive 
the turbulence and amplify the magnetic field. The j cr x B cannot extract energy 
from the CR because magnetic field only deflects particles and cannot change 
their energy. The solution here comes from the second order electric field CE2 = 
— ui xBi. E2 is anti-parallel to j cr o and reduces the CR energy which is transferred 
to the second order magnetic energy density TS 1 /8tt and kinetic energy density 
puf/2. 



4.4 The non-resonant instability in 3 dimensions 

In the case of a parallel shock, the zeroth order CR current is parallel to the zeroth 
order magnetic field and the above theory for k aligned with Bo can be applied 
directly. A more general theory is needed for an oblique shock where k and Bo 
are not paralle l. The disp ersion relation for general orientations of k, j and Bq 



isp< 
35) 



was derived by Bell (2005) for wavelengths shorter than the CR Larmor radius, ie 



12 K. M. Schure et al. 

kr g > 1: 

{ 7 2 + (k • hfk 2 v 2 A ) { 7 4 + l 2 k 2 (v 2 A + c 2 ) + (k • h) 2 k\ 2 A c 2 } (27) 

= 70 {l 2 + (J ' k) 2 fc 2 c 2 + [(k- b) 2 + (k- j) 2 - 2(k- j)(b- j)(k- b)] k\\] 

where v A = Bo/(4npo) 1 ' 2 is the Alfven speed, c s = (dP/dp) 1 ' 2 is the sound speed, 
7o = (k ■ Bo) 2 j 2 /p 2 c 2 , and a hat denotes a unit vector: k = k/|k|, b = Bo/|Bo|, 
j = jcr/|jcr|- The terms involving kv A are important at the short wavelength limit 
when magnetic field tension is important as in the case of k aligned with Bo . The 
terms involving kc s represent additional short wavelength compressibility effects 
that are not present when k, Bo and j CT are all parallel. The terms in kv A and 
kc s are important only at short wavelengths. At longer wavelengths they can be 
neglected and the growth rate simplifies to 



(k-Bo)j c 

PQC 



1/2 

(28) 



The instability grows most rapidly for wavenumbers parallel to the magnetic field 
but the growth rate is independent of the mutual orientation of the magnetic field 
Bo and the CR current j cr . 

The insensitivity of the growth rate to the angle between magnetic field and CR 
current implies that the instability is present for perpendicular as well as parallel 
shocks. In fact, the growth rate is faster for perpendicular shocks because the CR 
current is larger than that upstream of parallel shocks. In the precursor of a parallel 
shock the CR drift at the shock velocity relative to the thermal plasma to give a 
current density parallel to the shock normal of j C r = n cr eu s where n cr is the CR 
number density. In contrast, for a perpendicular shock the CR current density is 
larger in the direction perpendicular to the shock normal. The CR current density 
at a perpendicular shock can be calculated from the first moment of the Vlasov- 
Fokker-Planck (VFP) equation in the diffusive limit in which the one-dimension 
CR distribution function takes the form /(p, z,t) = /o(jp|, z, t)+fi(jp|, z, t)-(p/|p|). 

§ +c f>+^Ef>+^xf 1 = -„f 1 (29) 

at oz p op p 

where v represents angular scattering by small scale fluctuations. For a mono- 
energetic CR distribution at a perpendicular shock, the CR current density in 
the precursor can be separated into a component ju normal to the shock and a 
component jj_ that is perpendicular to the magnetic field and the shock normal: 

1 VOJg r g . 1 ULUg A . , 

•?ll = q^ - — ~-j-n cr ec J± = -^- — ncrec (30) 

11 6 oJg + v A L 6 tog + u z L 

where n cr is the CR number density, L is the scalelength of the precursor, A = c/v 
is the CR mean free path, and r g = c/uj g is the CR Larmor radius. 

The components ju and j± of the CR current density correspond to the CR 
drift velocities in the local upstream fluid rest frame. In steady state, ju = n cr eus 
and j± = (A/rg)n C reu s . The large scale field is relatively unimportant when the 
mean free path is comparable with the Larmor radius A ~ r g (Bohm diffusion) 
since the distinction between parallel and perpendicular shocks is then relatively 
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minor, and the scaleheight in both cases is L ~ (c/u s )r g . Bohm diffusion corre- 
sponds to the smallest possible mean free path. More usually, the mean free path 
can be expected to be greater than the Larmor radius and the cases of parallel and 
perpendicular shocks become quite different. From the above equations, the pre- 
cursor scaleheight ahead of a perpendicular shock is reduced to L ~ (r g /A)(c/u s )r g , 
and j± exceeds j\ i . j± results from the non-cancellation of gyratory currents in the 
CR density gradient in the precursor. Since jj_ is greater than j\ i , the instability is 
driven more rapidly at a perpendicular shock than at a parallel shock. However, 
the time during which the instability can grow is reduced because the scaleheight 
L is smaller at a perpendicular shock. The increased growth rate and the reduced 
time for growth cancel out and the number of linear e-foldings is the same for both 
parallel and perpendicular shocks. Hence the non-resonant instability is equally 
effective for all shocks whether they are perpendicular, parallel or oblique. From 
Equation 1281 in any of these cases fastest linear growth occurs for wavenumbers 
parallel to the large scale magnetic field, independent of its orientation to the 
shock normal. Unstable growth and magnetic field amplification at perpendicu- 
lar shocks has been demonstrat ed numerically in particle-in-cell simulations by 
iRiauelme and Suitkovskvl ( 201011 . 



4.5 Non-linear magnetic field amplification 

A linear instability takes a small perturbation 5B on the magnetic field and ampli- 
fies it until it becomes comparable with the zeroth order field Bo- When <5B ~ Bo, 
the linear assumption that second order terms in <5B can be neglected becomes 
untenable. A crucial question is whether the non-linear terms cause the instability 
to saturate and stop growing, or whether the magnetic field grows further to a 
magnitude much greater than Bo- In the context of CR acceleration, the field has 
to continue growing beyond SB /Bo ~ 1 if diffusive shock acceleration is to explain 
the presence of PeV CR in the Galaxy. Furthermore, saturation at SB / Bo ~ 1 is 
insufficient to explain the large magnetic fields inferred from X-ray observations 
of synchrotron emission at an SNR shock. Therefore it is crucial to ascertain not 
only that linear growth is sufficiently fast but also that the instability continues to 
grow non-linearly beyond SB/ Bo ~ 1 to generate fields exceeding 100 fiG at SNR 
shocks. 

Fortunately the non-resonant instability has the unusual property of continuing 
rapid growth into the non-linear regime. Remarkably, in the restricted geometry 
of a monochromatic circularly polarised wave with wavenumber k, zeroth order 
uniform magnetic field Bo and uniform CR current j cr all parallel, the linear equa- 
tions remain valid into the non-linear regime and the instability continues to grow 
exponentially to arbitrary amplitude at the linear growth rate. In this special case, 
slabs of plasma with frozen in magnetic field continue to be accelerated in direc- 
tions perpendicular to k, Bo and j C r- In practice other modes with different k 
also grow and these interfere to slow the growth. For example, an exponentially 
expanding spiral field in one part of the plasma is likely to collide with an ex- 
panding spiral field seeded in a different part of the plasma. The spirals cannot in 
general pass through each other and their growth is limited. However, the pres- 
ence of exponentially growi ng non-linear m o des is a strong hint that growth to 
large amplitude is possible. iLucek and Belli ( 20001) showed numerically that the 
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Fig. 2 Comparison of structures of density an d magnitude of magnetic field. 2D slice through 
a 3D simulation. Reproduced from [Belli 1120051) . 



magnitude of the magne tic fie l d can increase by at least an order of magnitude. 
3D MHD simulations by iBelll ( 20041) showed similar or larger growth before the 
calculation was terminated whe n magnetic structures expanded to the size of the 
periodic computational box. In IBelll ( 20041) the instability grew exponentially at 
the expected rate until SB /Bo ~ 1, whereafter it continued to grow but more 
slowly. 

The linear eigenmodes of the instability consist of spirals of magnetic with 
a preferred helicity. The evolutio n of this ba sic configuration into the non-linear 
regime can be seen in figure 4 of IBelll ( 20041) . Initially small spirals and loops of 
magnetic field grow non-linearly in radius. Collisions between neighbouring spirals 
produce walls of strong magnetic field surrounding cavities of very weak magnetic 
field. The field is far from uniform and does not conform to conventional pictures 
of randomly phased Fourier modes in fc-space. Because the magnetic field is frozen 
in to the background plasma, the density and magnetic field have closely corre- 
lated structures of walls and cav ities a s show n in Figure (2j which reproduces two 
frames from figures 2 and 3 of Bell (2005). Th e same wa ll-cavity struture has 
been found by Reville et al. (2008); Zirakashvili et al. (2008) in MHD simulations, 
and bv lRiauelme and Spitkovskvl (J2009I) and lOhira et al.l (J2009) in particle-in-cell 
simulations. 

The growth of large scale structures resulting from the expansion of the small 
scale structures provides a natural way of producing structures on the scale of a 
CR Larmor radius. These are especially important since CR are most effectively 
scattered by fields on this scale. Fields on smaller scales may explain the amplified 
fields observed at SNR shocks, but they cannot by the mselves provide the s trong 
CR scattering needed to accelerate CR to PeV energies. iReville et al.l ( 2008J) have 
modelled CR transport in non-linear CR-driven magnetic fields calculated with 
a 3D MHD code. They show that the amplified field does inhibit CR transport 
and reduces diffusion to less than Bohm diffusion in the initial magnetic field. 
The generation of magnetic structures on large scales is an active field of research 
(see section [5]), and other sources of turbulence contribute to the overall magnetic 
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field structure in SNR (Giacalone and Jokioii 2007; Zirakas hvili and Ptus kin 2008; 



Beresnvak et al.ll2009t llnoue et al.E o09: Schure et al. 2009) 



It is clear from analyses and simulations that magnetic field amplification con- 
tinues far into the non-linear regime. In some circumstances amplification may be 
limited by the time for which the j cr x B driving force operates. In the case of a 
shock precursor this limiting time is the time L/u s it takes for the precursor of 
scalelength L to be overtaken by the shock. However, two conditions have emerged 
as requirements for unstable growth. One is that the scale size fc _1 of magnetic 
structures should not exceed the CR Larmor radius, otherwise the CR follow the 
field lines and j cr x B becomes small since j cr is parallel to B. The second is that 
the magnetic tension B x (V x B)c/4-7r ~ kB 2 c/Aii should not exceed the j cr x B 
driving force. These conditions reduce to k > eB/pc and k < Atrjcr/Bc respectively. 
For both conditions to be satisfied simultaneously we need B 2 /An < pj cr /e. Using 
the expression for j cr derived in the discussion between equations [25] and [26] in 
section POl this reduces to an estimate for the saturated magnetic energy density 
producable in the non-linear phase of the non-resonant instability. 

**)«*" C HP2/P1) 

This estimate is a goo d match to observations of the magnetic field in SNR 



(Vink 2008; Bell 2009), but other explanations are possible for the same data 



(jMalkovet al.l l2011). 



5 Magnetic field amplification: long length scales 

As shown in the previous section, the non-resonant Bell instability can act effi- 
ciently to amplify magnetic fields on scales smaller than the gyroradius. However, 
in order to accelerate cosmic rays to higher energies, the magnetic field should 
also be amplified on scales beyond the cosmic ray gyroradius, which is the 'long- 
wavelength regime' discussed in this section. In 2011 various papers have been pub- 
l ished on possible instabilities that act to amplify the magnetic field on these scales 
(|Bvkov et al.ll2011bl: ISchure and Belll2011al: iReville and BdJl2012h . on which we 



will focus in this section. We will also briefly discuss other long-wavelength insta- 
bilities (firehose and acousti c) , but for a more extensive review refer to the original 
pape r s or reviews (see e. g. iBlandford and Eichlerlll987l : iMalkov and O'C Drurvl 
l200lt iBvkov et alJl2011ah . 



5.1 Current-driven stress-tensor instability 

The small-scale instability is essentially a fluid instability. The cosmic ray current 
is regarded to be of such scales that it is unperturbed by the non-resonant growth 
of small-scale magnetic fields. When looking at growth of magnetic field on scales 
larger than the gyro-radius, this assumption no longer holds since the CR tra- 
jectories follow the field lines and perturbations to the current in perpendicular 
directions should be taken into account. 

One way to determine this effect on the stability of the plasma is by including 
higher order anisotropic terms in the kinetic equation. The feedback between the 
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cosmic ray particles and the magnetic field in the plasma proceeds through forces 
acting on the momentum equation. The j x B forces arise when components of the 
cosmic ray current, or rather the induced return current in the plasma, are perpen- 
dicular to local components of the magnetic field. Gradients in the perpendicular 
current require taking into account higher-order components of the distribution 
function, at least the stress tensor, effectively representing gradients in the current. 
The distribution function of relativistic particles can be written as: 



d t { + v ■ Vf + — ( E + 12Lf±) . v„f = V„ • (D • V„f), (32) 

m \ c J 

where f is the particle distribution in phase space, v the particle velocity, and 

2 

D = ^^ (I — nn) the diffusion tensor with v the collision frequency, I the iden- 
tity matrix, and n the unit vector in the direction of the corresponding tensor 
component. 

The collisions, represented in the parameter v, are not actual collisions in 
such tenuous systems, but effectively act in the same way. What we mean with 
collisions is the cumulative effect of small-angle scattering as a result of the Lorentz 
force of the perturbed current and magnetic field. When they are fluctuating on 
the same scale, the Lorentz force is effectively deviating the path of the cosmic 
rays. Multiple of these small-angle scatterings result in isotropisation. The length 
and time scales on which this occurs is represented in the parameter v, which 
represents the scattering frequency in terms of the particle velocity over the mean 
free path. In Bohm diffusion this is taken to be of the order of the gyrofrequency. 
The short-wavelength non-resonant instability acts to amplify the field on scales 
that can efficiently deflect the low-energy cosmic rays. In effect, the parameter v 
thus describes the momentum exchange of the small-scale instability and can be 
used to determine its influence on the long range. 

The distribution of cosmic rays is dominated by the isotropic component: /o, 
but also contains an anisotropic part, to increasing order f = /o + f l • v/w + f 2 • 
vv/w + . . .. fi can be viewed as the directional component of the cosmic ray 
distribution, or as the gradient of /o, and acts like a current. f 2 is the pressure 
tensor, of which the isotropic part of the diagonal is normally included in the /o 
term. Anisotropy in the diagonal can be responsible for the firehose instability. 
Off-diagonal terms embody the stress-tensor and reflect gradients in the current. 
Each higher order is a factor of u/c smaller than the previous order, where u is 
the drift velocity. Evaluation of the transport equation to zeroth order, being the 
isotropic part, and first and second order anisotropies, gives the following system 
of equations, where we ignore any contribution of higher (3rd) order: 

9t/o = -|v-fi. (33) 

ftfi + cV/o + — (B x fi) + i/fi + |cV • f 2 = 0. (34) 

mc 5 



c I Vfi - \v ■ fil 2 ) + — (B x f 2 ) + 3i/f 2 

3 / mc 



= 0. (35) 

2 



Here I 2 is the second order unity tensor, and [. . .] 2 indicates a summation of the 
permutations for ijk in two ways divided by 2, such that we get a symmetric 
tensor with components that satisfy /^ = fji . In principle also higher order terms 
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can be used in the evaluation, but it t urns out these have no significant effect on 
the instability ( Schure and Bellll2011a|) . This system of equations can be closed in 



combination with the MHD equations: 

3(B = Vx(ux B) (36) 

pd t u = ith X B - VP - V ■ n + V Pcr (yL cr - u), (37) 

where p cr is the mass density of the relativistic particles, and u cr — u is the drift 
speed of the cosmic rays relative to the plasma, which to zeroth order is equal to 
u s . Since j t ^ = — j C r + c/(47r)V x B the momentum equation fEq. l37|) has to satisfy: 

pd tU = - jcr X B + — (V xB)xB-VP-V-i7 + U P cr(u cr - u). (38) 

C 4-7T 

In systems relevant for efficient diffusive shock acceleration, the upstream 
plasma can be considered cold and isotropic, such that VP = and V • H = 0. 
Additionally, the second term on the r.h.s. is much smaller than the other terms, 
and only contributes at very short wavelengths where the magnetic tension is suf- 
ficient to quench the instability. If we furthermore use that j C r = n cr eu s = efic/3, 
divide both sides by p, and write n = n cr /ni the ratio between cosmic ray- and 
background nucleons, the momentum equation can be expressed in terms of fi as 
follows: 

9tu = -^-(fixB) + |«mfi. (39) 

For a linear analysis it suffices to look at the first order perturbation, for which 
the above can be expressed, using Eq. 1341 as: 

9 t u (1) = -y (Vi ( i) + cV/ (i) + §cV • f 2(1) J , (40) 

where the subscripts between brackets indicate unperturbed (0) or perturbed (1) 
variables. 

Both of the two above equations are instances of the momentum conservation, 
viewed either through the forces j x B and frictional force pvu, or through the 
pressure gradient and divergence. Feeding this into the induction equation, we can 
express the perturbed magnetic field in terms of fi and f 2 : 

a t 2 B (1) = 9 t (B (0) • V)u (1) = -y (B (0) • V) (Wi(i) + cV/ 0(1) + |cV • f 2(1) Vl) 

where we assumed a homogeneous background magnetic field, incompressibility, 
and \ii \ = 0. For the rest of this section we consider a parallel shock in the z- 
direction, such that B( ) = P(o)^i an d we consider modes parallel to the original 

field, such that k ■ B( ) = 0. 

Equations I34IHT1 can then be combined fsee lSchure and Bel]||2011aj) to arrive 
at the dispersion relation: 
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Fig. 3 Growth rate of the tensor-mediated instability, from Schure and Bell (2011a). The 
solid (dashed) lines indicate the right- (left-) hand mode. The different colours give the growth 
rates for different values of the effective collisionality of the cosmic rays, which aids in coupling 
the cosmic ray momentum equation to the momentum equation of the plasma. 



where uj is the complex frequency, k the wavenumber, c speed of light, v the effec- 
tive scattering frequency, u) g the gyrofrequency, and fl = \/kjoB(n\/(pc) contains 

the driving (return) current jo and is the growth rate of the non-resonant Bell 
instability. The upper signs correspond to the left-handed polarisation (which is 
the polarisation of a gyrating proton), and the lower signs to the right-hand polar- 
isation. This dispersion relation is valid in the linear regime as long as the Alfvenic 
stress due to (V x B) x B is small (for the parameters used to plot Figure [3] beyond 
kr g « 1000). Also, since w is always small compared to fc 2 c 2 /(5(3^=F »w s )), it can 
in practise be ignored on the right hand side of Equation 1421 

The terms including factors of k 2 result from the inclusion of the stress tensor. 
Effects from the small-scale non-resonant instability are included through u, the 
effective scattering frequency. Scattering on small scales is expected to arise earlier 
than on long scales, since the growth rate increases rapidly for kr g ^> 1. Bohm 
diffusion is the regime where the effective scattering frequency is of the same 
order as the gyrofrequency, i.e. u w u g . We plot the growth rate as a function of 
wavenumber for different values of v/u g in Figure[3l for both the left-hand (dashed) 
and the right-h and polarisation (for the full plots of the dispersion relation see the 
original paper: ISchure and Bellll2011al) . 
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This method recovers the iBelll ( 20041) growth rate for the right-hand polarisa- 



tion in the regime where kc^> uj g . The resonant instability is only approximately 
captured in the mono-energetic approach (as can be seen from the peak around 
kr g ~ 2 in Figure [3]) , which is why the familiar fc-dependent growth rate around 
kr g = 1 is not present in this representation. Care should be taken when including 
the resonant instability to do the momentum integration with an appropriate up- 
per limit for p, so as to not overestimate the growth rates on the long-wavelength 
end. It can be seen that both modes are unstable. Which of the polarisations 
dominates depends on the ration v/u)g\ for v /uj g < 1/vo the left-hand mode dom- 
inates, and vice versa for higher collisionality. When the collisionality is zero, 
v = 0, only the left-hand mode is unstable and purely growing. The current-driven 
long- wavelength instability depends on the mediation of the short-scale instability 
through the stress-tensor and the 'collisionality'. 

Self-consistency between the equations for the cosmic rays and for the fluid are 
crucial in deriving the dispersion relation. In the following analysis we show that 
omission of the force due to friction in the momentum equation, as can often be 
done in other circumstances, would result in a completely different dispersion rela- 
tion that shows a much more rapid growth and only declines at long wavelengths 
as yk. This is a result of an unbalanced frictional force, that results in the vi\ 
term remaining in the momentum equation (Eq. I40[) as a consequence of there not 
being a similar term in the momentum equation for the fluid. This results in the 
following erroneous dispersion relation: 



.2 i n2 



w 



7 Z A \ / 7 2 2 

k c \ . I . . k c 



0{6U^lOJg) J \ 5(3/V=FlWg) ' 



The difference is only an additional v in the nominator. For large k, the additional 
v can be ignored such that on the short scales the result does not change. However, 
on the long- wavelength end the additional v would change the result. To lowest 
order in k (thus ignoring the k contribution from f2), and in the limit w <JC ui g , 
the growth rate would change to: 



N 



2(^ 2 +a;2) 



(44) 



In the limit where v <g; ui g , both modes would have similar growth rates, which 
would read: 



It should be stressed this is not a physical solution 
5.2 Ponderomotive instability 



^=J^T^ = \I^T ("«"»)■ ( 4 5) 



iBvkov et al.l ( 2011bl) presented a long- wavelength instability that results from av- 



eraging the kinetic equation for the relativistic particles, the equations of the bulk 
plasma motions and the induction equation over the ensemble of the short scale 
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fluctuations produced by CR instabilities in the collisionless regime e.g. by the fast 
Bell instability. To derive the growth rates of the modes in the long-wavelength 
regime kA < 1, with A = r g /(y/ui g ) the mean free path, the dispersion relation as 
in Eq. 1611 derived from the collisionless kinetic equation approach is in this case 
not appropriate. 

In the presence of the short-scale fluctuations, the momentum exchange be- 
tween the CRs and the flow in the hydrodynamic regime, results in a ponderomo- 
tive force that d epends on the CR c urrent in the mean-field momentum equation 



of bulk plasma ( Bvkov et al.ll2011br) . As a result, there exist transverse growing 



modes with wavevectors along the initial magnetic field with growth rates that are 
proportional to the turbulent coefficients determined by the short scale fluctuation: 



2B pcuj g ' [W) 



where \/{b 2 } is the magnitude of the short-scale amplified magnetic field and which 
holds for both polarisations. The magnetic field amplification in that regime only 

— 1/2 

weakly depen ds on the shoc k velocity (-yr oc u s a s follows from Eq. (48) 



weakly depends on tne snocK velocity (jr oc u s ' as lollows Irom fcq. (48J in 
iBvkov et al.l ( 2011bl) . see also ISchure and Belli ( 2011b|) ]. that is important for the 



evolution of the maximal energy of CRs accelerated by DSA. In the intermediate 
regime, v jio g < kr g < 1, the growth rate can be approximated as: 



7 = 47rfcv / <v 5 ), (47) 



where \/(v 2 } is the amplitude of the short-scale turbulent bulk velocity. 

The ponderomotive instability is a multi-layered phenomenon and the under- 
lying physics is not immediately clear. In the intermediate regime, the growth rate 
(Equation I47|) is independent of the CR current which suggests that it is not di- 
rectly driven by CR streaming. Instead, the growth time is equal to the time taken 
to cross a distance 1/fc at the characteristic turbulence velocity. This suggests that 
the magnetic field in this regime grows as a result of field-line stretching by al- 
ready existing turbulent motions. The long wavelength regime (Equation I46[) is 
complicated because it includes magnetic field on three different scales: i) Bq on a 
scale comparable with or greater than the wavelength, it) the magnetic field \/{h 2 } 
associated with the turbulence, and Hi) the small scale magnetic field causing the 
scattering represented by the collision frequency v. The challenge is to see why all 
three fields are important for the instability. Is also important to check that the 
j x B momentum transfer between CR and the fluid is treated self-consistently in 
each case since errors in total momentum conservation can lead to an incorrect 
growth rate as shown in section 5.2. In the next section we review the filamenta- 
tion instability which is also due to the presence of turbulent magnetic field on 
scales smaller than kr g ~ 1. A further challenge is to ascertain whether there is 
any overlap in underlying physics between the ponderomotive and filamentation 
instabilities. 



5.3 Filamentation instability 

The expansion of the loops generated by the non-resonant instability on small 
scales can give rise to a further filamentation instability. Because the cosmic rays 
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are focussed into filaments, the cosmic rays current locally increases. As a result 
magnetic fields around the loops further grow in strength, which again aids to focus 
the cosmic rays and incre ase the current. This w as recently derived analytically 
and shown numerically bv lReville and Belli ( 20121) . The growth rate turns out to 



be independent of wavelength. Again, the Vlasov equation is used to determine 
the distribution of cosmic rays. The local electric field can be expressed in terms 
of the vector potential such that: 

E = u s VA||/c, (48) 

with An the magnitude of the vector potential parallel to the shock normal. Using 
this, the distribution function can be written as: 

% + cH ■ V/ + eV(« fl A« /c) • ^ = 0. (49) 

at p " dp 

Since it can be safely assumed that df /dp < 0, the cosmic-ray number density is 
locally larger when An is positive and can be written as a number of position: 

n C r = n -\ 2 — / 87r P/o d P, (50) 

where no = J Anp fodp and /o the isotropic part of the distribution function. 
Using further that j C r = n cr eu s and the MHD equations, the evolution of the 
filamentation can be expressed in terms of the cosmic ray current as follows (see 
iReville and Bellll2012h : 

d 2 j CI e 2 n Q u 2 .B ± . _..,_,. , T7\ 5 Jcr , k1 n 

^^ = ^r^^ Jcr + ([(u - v)u] ' vbcr_(u - v) ^r- (51) 

The first term on the right-hand side is independent of the wavenumber and ob- 
viously dominates the other terms on long scales. Its growth rate is: 



7 = a/^ = r-) 2 (%y /2 *± , (52) 

c y pc p x \ c J \puij piy/ln(p 2 /pi) 

for a cosmic ray spectrum with a power law slope q = 4 with a given minimum 
(pi)and maximum (p 2 ) momentum cut-off. The requirement on the cosmic rays 
driving the instability is that they are not trapped within the cavities, such that 
p\c ~^> eA«Us/c, with An the vector potential. This condition was also assumed 
when deriving the instability. The growth rate is linearly dependent on the value 
of the amplified small-scale magnetic field and decreases upstream of the shock 
when pi increases and (7 C r decreases. 

When compared to the non-resonant growth rate on small scales, the growth 
rates are equal when 

*r. = ^<§>. (53) 

Non-linear simulations of the no n-resonant instability indicate that the am plified 
field can reach values of 30Bq ( Belli l2004t iRiquelme and Spitkovskvl 120091) . such 
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that the above condition can be satisfied for kr g > 1. Filling in numbers comparable 
to those used in the previous sections, the instability operates provided that: 

E - >~ (life) (jSlfc) _1 (lO^nTi) ^ ^ 

When compared to the growth rate of the long- wavelength instability, they become 
comparable when: 

Thus, for the longest wavelengths, the filamentation instability may dominate if 
the small-scale field is sufficiently amplified and as long as conditionl54lis satisfied. 



5.4 Firehose instability 

The short-wavelength non-resonant instability in Section 14.21 was driven by the 
current, i.e. the first order anisotropy in the cosmic ray distribution. Asymmetry 
in the second order anisotropy, the pressure tensor, can drive the firehose or mirror 
instability. It is distinctly different from the long-wavelength instability that was 
discussed in Section [5TT] that one was driven by the current, and the stress-tensor 
(off-diagonal terms of the pressure tensor) and small-scale collisions mediate this 
driving source to cause instability on long length-scales and the opposite (left- 
hand) polarisation. 

Again, we will have to consider the distribution function up to second order 
anisotropy, which we can write in the form 



K \ N 

/(P) = 4^ 



1 + 3^ + 1(3^ 



(56) 



where 9 - particle pitch-angle, /i = cos 8 , S(p) - is the magnitude of the second har- 
monic anisotropy, which is normally of the order of u 2 /c 2 . Indeed, it is the second 
harmonic anisotropy that constitutes the source of the CR-firehose instability on 
the magnetic field amplification. It is instructive to summarize the growth rates for 
magnetic instabilities that the quasi-linear theory predicts for weakly anisotropic 
CR distributions of the above form. Following the standard linear a nalysis of the 
kinetic equation in the intermediate regime r g j A < xo < 1 (see e.g. iBvkov et al.l 
12011 q) . with A again being the mean free path, one may get the following disper- 
sion relation if collisions can be neglected: 



vlk 2 



i*S{c*-u + £* 



(57) 



f P2 2 

A 0A (xi,x 2 ) = / ao,i (p) N (p) p dp (58) 



-^=!£w^ (59) 



MFA and DSA 23 

<J1 ^ = 4,/_ 1 1 T „ ^ (60) 

where fc = S— ^ as = kr g(p) > x i = kr g (pi) j »2 = fcr g (p 2 )- The signs ± 

C -DO 

correspond to the two opposite circularly polarized modes under investigation. 
The second term on the right hand side is the instability that is driven by the 
cosmic ray current, and has a non-resonant form (|Bellll2004l , see previ ous section) 



when kr a is large, and a resonant form around kr a of the order (1) (e.g . Achterbere. 



1981 



; Zweibel 2003; P elletier et aTll2006l iMarcowith et al.ll2006t lAmato and B lasi 



2009, and the references therein). The last term in the r.h.s. of Eq. 1571 represents 
the CR firehose instability. In the long-wavelength regime x m <C 1 a simplified 
form of Eq. [5jJ can be derived 



vlk* l T 



koTgO 



Sc/u s \n(p 2 /pi) 
X2 (l-^M)" 1 ) 



(61) 



The growth rate of the firehose instability due to the CR pressure anisotropy is 
found in the last term of Eq. [ST] It requires that the parallel pressure exceeds 
the perpendicular press ure such that Pn > P± + B 2 /(An) (for more details see 



tne perpendicula r press ui 
e.g. iBvkov et al1l2011ah . 



5.5 Instabilities driven by the cosmic ray pressure gradient in the shock precursor 

On scales large compared to the gyro-radius (i.e. scales where the driving particles 
are strongly magnetised) the action of the cosmic rays can be simplified to a bulk 
force on the background plasma which is just given by the gradient of the cosmic 
ray pressure, 

Per = J ^p 2 /(p)dp (62) 

This follows naturally from considerations of momentum balance, the above pres- 
sure integral being just the flux of momentum associated with the cosmic rays. 

The interesting thing about this bulk force is that it is not related in any 
simple way to the local density either of mass or of scattering centres (this is a 
consequence of the collective nature of the electro-magnetic interactions in the 
plasma). A gravitational field, by contrast, would produce a force that is always 
strictly proportional to the local mass density (this is just Einstein's equivalence 
principle) . Similarly a flux of particles interacting by two-body scatterings would 
produce a force that is proportional to the density of scattering centres, radiation 
pressure resulting from Thompson scattering on electrons being an example. No 
such simple relation holds for the cosmic ray scattering, which is a complicated 
function of the power-spectrum of structure in the magnetic field on scales compa- 
rable to the particle gyro-radius and, to lowest order, can be calculated using the 
methods of quasi-linear theory. For our purposes it is enough to just assume that 
there is some effective diffusion coefficient n(p) and that on the scales of interest 
to us the cosmic ray transport can be represented as a simple diffusion process 
with a flux proportional to the local gradient, 

^ + «-V/=iv-« P | + V.( K V/) (63) 
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where k is the diffusion tensor (often approximated as a scalar diffusion coefficient) . 
A given element of background plasma, of local density p experiences an ac- 
celeration 

|£ = -±VP cr (64) 

at p 

and thus, even if the cosmic ray pressure is very uniform, local small-scale vari- 
ations in density induce acceleration fluctuations which in turn lead to velocity 
fluctuations which can feed back into density fluctuations. The usual approach to 
studying instabiities, where one assumes a uniform steady background and then 
does a Fourier analysis of the modes, fails in this case because the non-stationary 
and non-uniform nature of the shock precursor region is the ultimate source of 
free energy driving the instability. Indeed the very question of how to define an 
instability in such a system is an interesting one with no obvious answer. 

In iDrurv and Falld ( 1986|) a solution to these problems was developed based 



on a two-scale expansion of the governing equations. In this approach one looks 
at small wave-length high-frequency modes propagating on a smoothly varying 
background. In the absence of the cosmic ray pressure the basic modes are then 
just sound waves and the wave amplitude satisfies a conservation equation for the 
wave action (this follows from Noether's theorem because in the high-frequency 
limit the precise phase of the wave becomes unimportant, and there is thus an 
asymptotic symmetry related to the arbitrariness of the phase angle). In the case 
of one spatial dimension this is 

^ + §- z [(u±c s )A}=0 (65) 

where the wave action A is the acoustic wave energy density divided by the local 
co-moving frequency, c s is the sound speed and the sign of ± corresponds to left- 
and right-travelling modes. 

Inclu ding cosmic ray effects i n the two-fluid approximation it is then possible to 
show fsee lDrurv and F alle 1986) that the wave action equation acquires a non-zero 
right hand side, 



- + - [{u±Cs) A} = A 



JcrPcr , dPcr/dz f dlnK 



pn '" c s p \ dlrip 



(66) 



where y cr is an effective adiabatic index for the cosmic rays and dlnn/dlnp is 
the extent to which fluctuations in density induce fluctuations in the diffusion 
coefficient. 

The first term is a linear dam ping term related to the cosmic ray diffusion, as 
derived earlier bv IPtuskinl ( 19811) . More interesting for our purpose is the second 



term which is a potentially de-stabilising term related to the cosmic ray pressure 
gradient. By formulating the problem in this way it is possible to clearly sepa- 
rate out the conservative effects of the changing background, encapsulated in the 
conservation of wave action, from the non-conservative effects of the cosmic ray 
pressure. This gives a precise instability criterion, 

1 dPcr ( \ , d\nn\ ■JcrPcr , c „\ 

1 + 7TTT7 > —— ( 67 ) 



pc s dz \ dlnpj pn 
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A gradient in the cosmic ray pressure can arise when the cosmic rays make up a 
significant fraction of the total pressure, and is of interest in the shock precursor 
region of an efficiently accelerating shock wave. 

If we introduce a length-scale for the cosmic ray pressure, 



^ = P ~f < M > 



the condition for instability can be written as 



L< _*_r 1+ o£*\ (69) 

jcrCs \ amp) 



Noting that in a shock precursor L ~ «/w s where u s is the shock velocity, it is 
clear that the shock precursor region will be generically unstable at high shock 
Mach numbers unless <91nK/<91np is identically — 1. This is confirmed by numer- 
ical simulations of modified shocks in the two-fluid approximation which exhibit 
instabilities unless the product pK is artificially kept constant. Indeed it was this 
phenomenon in the early calculations of Dorfi which led to the discovery of the in- 
stability. It should be noted that the instability only operates as a fluid element is 
advocated through the precursor, so the maximum growth is limited to an amount 
of order exp(A/) which can however be very large for high Mach numbers. 

In addition to the acoustic modes there are also non-propagating entropy 
modes. The entropy of each fluid element is conserved and therefore the entropy 
modes can not grow, at least until secondary shocks form. The modes do couple 
to acoustic modes which can then be amplified. Such entropy modes are perhaps 
better thought of as dense clumps in pressure equilibrium with their surroundings 
and the differential acceleration forces resulting from the cosmic ray pressure will 
set these in motion relative to their less dense surroundings. This motion will then 
be transmitted to the surroundings in the form of acoustic waves which, when 
propagating parallel to the shock normal, will be amplified unless the diffusion 
coefficient scales inversely with the density according to the above analysis. If we 
now consider the very special case of a sinusoidal density perturbation with wave 
vector perpendicular to the shock normal, it is clear that it too will induce motions 
unless the diffusion is constant and independent of density. Thus in three dimen- 
sions it is impossible to fully stabilize a clumpy shock precursor. The condition for 
parallel stability implies transverse instability and vice- versa. 

It is also worth noting, as pointed out bv iGiacalone and Jokipiil ( 20071) . that 



even with no precursor effects a clumpy upstream medium will induce strong post- 
shock vorticity and down-stream magnetic field amplification. The above analysis 
indicates that similar processes can work upstream if there is a strong cosmic- 
ray precursor and thus produce magnetic field amplification by what is in essence 
a bulk hydrodynamic effect. The scales on which this takes place can be large 
compared to the particle gyro-radii and is determined by the characteristic length- 
scales of the initial density fluctuations. 
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6 Deviations from the "-2" power law index 

6.1 The spectral index at oblique shocks 

The relative motion of upstream and downstream scatterers imparts energy to CR 
as they bounce back and forth across a shock. By Lorentz transformation, the mean 
fractional energy gain is u s /c on each passage from upstream to downstream and 
back to upstream for a strong non-relativistic shock. As discussed in Section[2l the 
resulting CR spectrum is determined by the statistical distribution of the number 
of times a CR crosses the shock before escaping downstream. CR cross the shock 
at a rate n s c/4 where n s is the CR number density at the shock. This is balanced 
by the rate noo Ms /4 at which CR advect away downstream of the shock, where noo 
is the CR number density far downstream. Hence the fractional number of CR lost 
after each crossing is (rioc / n s ){u s / c) . In the limit of small shock velocity, the VFP 
equation dictates that noo = n s and the fractional number lost is u s /c. When 
combined with the fractional energy gain of u s /c, a power law CR distribution 
results with differential energy spectrum n(E)dE oc E~ 1 dE where 7 = 2. If, as we 
show below, n s can differ from noo, the spectral index is 

7=1 + noo/n s . (70) 

The above argument predicts a universal E~ spectrum for all high Mach num- 
ber non-relativistic shocks, whether they are perpendicular, parallel or oblique. 
However, the argument depends on the result that noo = n s . It has been clear for a 
while that this breaks down when the shock velocity is relativistic (e.g. [Achterberg et al.l 
l200ll and Spitkovsky 2011 these proceedings), but it has recently been shown (Bell 
et al 2011) that departures from this spectrum occur at shock velocities as low 
as ~ c/30 as found in young supernova remnants (SNR). The departure is par- 
ticularly strong for shocks that are nearly perpendicular and when the CR mean 
free path A is larger than the Larmor radius r g . As shown in the previous sec- 
tion, the precursor scaleheight at a perpendicular shock is L ~ (r g /A)(c/u s )r g . 
According to this formula, the upstream scaleheight can be very short causing a 
discontinuity in the CR density gradient across the shock. Kinetic theory does not 
allow discontinuities occur over distances less than a CR Larmor radius since CR 
gyration imposes a smoothing distance of a Larmor radius on the CR distribution 
function. If the precursor scaleheight is not much larger than the Larmor radius, 
the overall CR change in density across the shock takes place partly downstream 
as well as upstream. The result is that the CR density n s at the shock is less than 
the density far downstream noo , the spectral index 7 is greater than two, and the 
CR spectrum is steepened as indicated by equation 1701 

In contrast, solution of the VFP equation shows that the spectrum is flattened 
if the shock is oblique and more than 10—20° from perpendicular, depending on the 
shock velocity. Compression at the shock increases the perpendicular component 
of the magnetic field because it is frozen in to the background plasma, whereas 
the parallel component of the field is unchanged by the shock. Consequently the 
magnetic field increases in magitude and changes direction at an oblique shock. 
The change in field acts as a partial magnetic mirror which reflects back upstream 
some CR trying to cross the shock into the downstream plasma. The shock acts 
as a partial snowplough, pushing CR ahead of it. This produces a local excess in 
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the CR number density at the shock such that n s > noo , and the CR spectrum is 
flattened in accordance with equation 1701 giving 7 < 2. 

Clearly, the universal strong shock spectrum, 7 = 2, does not hold for young 
SNR shocks. This is is consistent with radio observations which ex hibit significantl y 
steepened spectra in very young SNR expanding at high velocity ( Bell et al.ll201ll) . 
The effect is not confined to high velocity shocks. The crucial parameter is the ratio 
of the shock velocity to the velocity of the accelerating particle, so the spectra of 
sub-relativistic particles accelerated by heliospheric shocks may be expected to 
show a related departure from 7 = 2. 



6.2 Non-linearity and time-dependence 



From observations only direct evidence of acceleration of electrons in supernova 
remnant blast wave exists through the observation of narrow synchrotron rims in 
X-ray and radio-synchrotron emission from a more extended region. In theory there 
is no reason why protons and heavier nuclei would not be accelerated through the 
same process. Gamma-ray observations are currently hinting towards observations 
of escaping nuclei interacting with molecular clouds, and perhaps also in situ pro- 
tons at SNR shock waves. Protons could potentially be much more abundant (esti- 
mates of a factor 1000 are not uncommon), due to easier injection and less radiation 
losses, and may reach higher energies. Efficient proton acceleration would therefore 
have important consequences on the shock structure and temperature. The non- 
linear effects due to cosmic rays constituting a signific ant fraction of the energy at 



the shock have been widely discussed and modelled ([E llison et al. 1995; Malkov 
1997tlBlasill2002l:jKang and Jonesll2005t [Amato and Blasill2006HVladimirov et al 



2008 : lKang et alil2009i IPatnaude et al.M2009l: ICaprioli et alJl20ld: iFerrand et al 



2010) . Especially towards the higher-energy end of the spectrum, the spectral index 



can significantly flatten due to the higher overall compression ratio that is probed 
by the more energetic cosmic rays. On the low-energy end spectral steepening can 
occur, although this is only important in l ow-Mac h number shocks unless almost 
all of the energy goes into cosmic rays fe.g. lVink et al.ll2010T ). It may be important 
to consider in simulations of cl usters of galaxies, in which low-Mach number shocks 
appe ar to accelerate particles ( Rvu et al.ll2003c IPfrommer et al.ll2006clVazza et al.l 
I2009T) . The transfer of energy to cosmic rays may further be seen in temperature 
deviations behind the shock, where part of the energy that w as supposed to heat 
the plasma has effect ively gone into a cosmic ray component ( Helder et al.ll2009t 
IPatnaudeet al.ll2009l) . 

Whether a shock can be an efficient accelerator, and what the resulting cosmic 
ray spectrum looks like, is also dependent on the environment the shock is running 
into. A core-collapse SNR will have a different evolution of the shock velocity from 
a SNR evolving in a homogeneous medium as may be the case for most Type 
la SNe. The time-dependent evolution will determine the cumulative spectrum 
( Schure et el1.I I2010T) . The environm ent may in some ca ses also affect the damp- 
ing through ion- neutral collisions ( Reville et al.| |200J). the de tectable emission 
through various energy-exchange processes ('e.g. lRavmond e t al. 2011) or through 
surrounding moldecular clouds enhancing the target density for pion creation from 
escaping cosmic ray protons (e.g. lOhira et al.ll201ll) 
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Apart from magnetic field amplification, progress has been booked on the the- 
ory of DSA that deviates from the ideal case. The powerlaw of the cosmic rays in 
reality is a complex addition of time-dependence, shock obliquity and shock speed, 
as has been discussed in Section 16.11 The theory on how cosmic ra ys may escap e 



as has been discussed m section lb. II 1 he theory on now cosmic ra ys ma y escape 
upstream is als o unde r active development, with recent papers bv iDrurvl ( 20111) ; 
lOhira and Iokal (|201lh 



7 Discussion and conclusion 

The importance of magnetic field amplification has been recognised since the early 
developments of the theory of DSA, and in the past decade significant progress 
has been made in this field. Driven by observations that strongly indicated am- 
plification of factors 10-100, various theories have been developoed to explain this 
intrinsically nonlinear process. A distinction between scales shorter and longer 
than the gyroradius of the driving cosmic rays is made. Most rapid amplification 
can be achieved on short scales, but amplification on longer scales is paramount 
in accelerating cosmic rays to the PeV energies that they are believed to gain in 
galactic sources. Currently, SNRs still seem the best candidate, but the process is 
independent of the type of shock wave and other sources may well contribute. The 
nonlinear behaviour of the various instabilities remains an area of active research 
and more work is required to satisfactorily decide in which regimes various insta- 
bilites operate and dominate. Also the saturation level needs to be determined. 
Observations that determine more precisely direction, degree of polarisation, on 
different length and time scales could aid in constraining the theory. The high- 
energy end of the emission spectrum can be better probed with the current ob- 
servatories that operate over a range of wavelengths in the gamma-ray regime. 
Most notably the low-energy end as probed by Fermi-LAT, as well as Cherenkov 
telescopes for the high-energy end, such as HESS. The next generation telescopes, 
such as CTA, will certainly be extremely useful, if the observed energy is pushed 
up high enough to really distinguish between electron- and proton- based emis- 
sion processes. High-resolution radio measurements of magnetic field strength, and 
polarimetry in the X-ray band, are other items on the wish list. 
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